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Following our paper [J. Math. Phys. 50 (2009) 123102], we systematieally earry out Lie symmetry 
analysis for the barotropic vorticity equation on the rotating sphere. All finite-dimensional sub- 
CsJ ■ algebras of the corresponding maximal Lie invariance algebra, which is infinite-dimensional, are 

classified. Appropriate subalgebras are then used to exhaustively determine Lie reductions of the 
equation under consideration. The relevance of the constructed exact solutions for the description 
of real-world physical processes is discussed. It is shown that the results of the above paper are 

■ directly related to the results of the recent letter by N. H. Ibragimov and R. N. Ibragimov [Phys. 
. Lett. A 375 (2011) 3858] in which Lie symmetries and some exact solutions of the nonlinear Euler 

^ : equations for an atmospheric layer in spherical geometry were determined. 

I ; 

1 Introduction 

One approach to reduce the complexity of the governing equations of geophysical fluid dynamics 

■ is to consider submodels thereof. The latter usually take into account several simplifications that 

■ make them adapted for distinct physical processes or that render them valid only on well-defined 
spatial or temporal scales. One of the most relevant simplifications which can be reasonably 
implemented in the atmospheric sciences is that of the two-dimensionality of large-scale flow. 

PsJ ■ This reduction in dimensionality greatly simplifies the structure of the dynamical part of the 

. governing equations, which can be consistently described by some forms of the two-dimensional 

Euler or Navier-Stokes equations [21 O El El [TH [TB] . Though they are still coupled systems of 

■ nonlinear partial differential equations, much is known by today about the long-term behavior 
of these systems. 

CN \ A form of the two-dimensional Euler or Navier-Stokes equations which is of immediate in- 

terest in large-scale atmospheric dynamics is that in spherical geometry. In the recent letter [8] 
the system of nonlinear incompressible Euler equations in a single thin atmospheric layer on the 
sphere rotating with constant angular velocity was considered, 

^ ■ u „ V cos 9 1 



Ut + Ux + vug + VU cot 6 H h : P\ = 0, 

sm a Hn sm 



U 2 /I " '^'-'^ 

vt + vx + VV0 -u cot e \-p0 = 0, 

sm a Kq 

Ux + {vsm9)g = 0. 

In the above system A and 6 denote the azimuthal and the polar angles, respectively. The former 
increases by going East, while the latter increases by going South. The associated velocity 
components are u and v, respectively, p denotes the pressure and Rq is the Rossby number, 
which is the ratio of inertial to Coriolis force, Rq oc The last equation of the above system 
is the incompressibility condition. Subscripts as usual denote the partial derivatives. 
In [8] it was assumed that 

u = F{9) + VU, v = vv, p = p + up, (2) 
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where is a (not necessarily small) parameter, F(9) is an arbitrary function of its argument 
and p is the mean pressure in the thin atmospheric layer studied. Variables with hats are to be 
interpreted as disturbance quantities. Physically, this ansatz corresponds to finite disturbances 
that are superimposed to a zonally averaged mean flow having only a latitudinal dependence. 
By introducing a stream function ■0 for the disturbance velocities, i.e. 

U = ljjg, V = T^i'X, 

smt) 

and taking the curl of the first two equations of the system ([1]), one arrives at 

C. + IM. - + C. - -i^ t.L.F + i fa = 0, (3) 

where the vorticity C is related to the stream function ■0 through the Laplacian in spherical 
geometry, 

C := {sin 94)0)0 + -\- V-AA, 
sm sm^ 6 

and Li denotes the Sturm-Liouville operator for the associated Legendre functions, 

T 1 d / d\ 1 

Li := — — -— sm( 



sin^d^V sin^e' 

This is the equation considered in [8]. 

Previously, in [Ij we have studied the Lie symmetries and exact solutions of the well-known 
barotropic vorticity equation on the rotating sphere, 

Ct + ^aCm - V'a.Ca + 2J^0A = 0, C := + - (4) 

which is also derived from the system of nonlinear Euler equations ([T]) , but without introducing 
the ansatz ([2]). In the above equation the variable /u is defined as fi = sin 93 with ip being the 
geographic latitude on the sphere. The mean radius of the Earth is assumed to be scaled to the 
unity. 

It is the purpose of the present letter to relate the above two independent studies. In Section [2] 
we explicitly transform Eqs. ([3]) and to each other. This allows us to derive the results on the 
Lie symmetries and reductions of Eq. ^ obtained in |8] from our previous work [3] . The related 
Lie reductions and exact solutions of Eq. Q are discussed in Section [H Additionally, in Section [3] 
we present the exhaustive classification of finite-dimensional subalgebras of the maximal Lie 
invariance algebra of Eq. (jH) . In Section \5\ we give some concluding general remarks on Lie 
reduction and illustrate them with the reductions carried out for the nonlinear Euler equations 
in a single thin atmospheric layer. 



2 Relation between the equations 

It is straightforward to show that Eq. Q and Eq. can be mapped to each other by point 
transformations. In fact, their difference merely comes from a different set of variables employed 
and the special ansatz ([2]) adopted in [8]. 

First we should like to note that in geophysical fluid dynamics it is common to use the latitude 
(p rather than the polar angle 9 as an independent variable. They are related by 99 = 7r/2 — 9, 
hence /i = sinip = cos0. The second difference between the two forms of the vorticity equation 
stems from the ansatz ^ assumed in [S]. This ansatz plays no essential role for the study of 
nonlinear Euler equations as it is explicitly assumed that 1^ is not a small parameter. Therefore it 
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is possible to avoid the separation of F{9) from the disturbance velocity v and the introduction 
of the parameter expansion in i/. As shown in Section 21 this ansatz even complicates the 
structure of exact solutions that can be derived. Solutions in the form of a separated mean flow 
superimposed by some disturbance can be generated by the exact solutions that can be derived 
from the usual form of Eq. In fact, setting 

= ^ (^V - y Fde^ (5) 

and using the triple of independent variables {t, A, fi) rather than (t, A, 9) leads to the cancelation 
of the terms involving F and thus to the classical barotropic vorticity equation in a rotating 
reference frame Q. Note that the above transformation is singular for = 0, a case which yields 
an exact solution of the incompressible Euler equations for F = 1/ sin 6 as discussed in [8j . 
A further simplification of Eq. (jj]) can be achieved upon using the transformation 

i = t, fj- = fJ-, A = A + Ot, ip = ip — Qfj,, (6) 

which was first found in [13] and then rediscovered in [3] using methods of group analysis. 
This transformation allows mapping the vorticity equation in a rotating reference frame to the 
corresponding equation in a reference frame at rest, i.e. to set = in Eq. 

3 Subalgebras of maximal Lie invar iance algebra 

Taking into account all the transformations discussed in Section [21 the maximal Lie invariance 
algebra presented in [5] is of course isomorphic to the maximal Lie invariance algebra of the 
vorticity equation ^ for = 0, which is generated by the vector fields [1] 



sin A 



V = tdt- ipd^, dt, Ji = dx, J2 = fi^==dx + vT^-^cos Xd, 



^^^^ (7) 

^3 = lJ'^==d\ - \/l - sinA(9^, Z{g) = g{t)d^, 

vi -/^ 

where the parameter g traverses the set of smooth functions of t. The algebra has the 
structure of so(3) ® (02 ^ {^ig)))i it is the direct sum of so(3) (with the basis elements J7i, 

and iTs) and the semidirect sum of the two-dimensional non-Abelian algebra 02 = {T^,dt) 
with the infinite-dimensional Abelian ideal {2(g)). Eq. @ admits also two independent discrete 
symmetries, (t,X,fi,ip) 1-^ {—t, —X, fi,tjj) and (t,X,fi,ip) 1-^ (t, X, —fj,, —ip). 

In order to systematically carry out Lie reduction, an optimal list of low-dimensional subalge- 
bras of the maximal Lie invariance algebra is required, which was also determined in [Ij. As 
subalgebras of greater dimensions may be useful, e.g. in the study of partially and differentially 
invariant solutions [21 [3 [HI [T2] and invariant parameterizations [31 \T0\ [T^ [T7] , here we classify 
all finite-dimensional subalgebras of . Up to equivalence generated by the internal automor- 
phisms of , a complete list of such subalgebras is exhausted by the following parameterized 
classes of subalgebras: 

1) 2) {Ji+Z{f)) ffi rig), 3) {Ji, J2, J3) © 

4) {d,+aj,) © i^-,^„ 5) {d„ J^+Z{Kt>^)) © /£^'-,^„ 6) (a„ Ji, J2, J,) © 

7) {V+aJ,) © I£^',,^„ 8) (P, J^+Z[nt-\^)) © i^^^, 9) (P, J,, J2, J3) © /S,,^,, 

10) {v+aji+zi^e), dt) © 11) {V, dt, ji+z{^e)) © i", 12) {v, dt, Ji, J2, Js) © 

where for convenience we use several notations for subalgebras contained in the ideal {Z{g)), 
r{-g) = {Z{g'),...,Z{gn), 
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= {Z{r^'\A'% Zirh\tp cos(i^r)), Z(r'^\tp sm(z^r)), l^^-:^^, Jlt;^)' 
I" = (2(l), 

Here ^ = (5^, . . . , (7*^) is an arbitrary n-tuple of linearly independent functions of t, n,n' G Nq, 

is an n-tuple of different real numbers, /i and v are n'-tuples of real 
numbers such that all pairs {iJ,^,^^) are different, r = ln|t|. In the fifth (resp. eighth) class we 
have k = ki with the value of i such that Aj = (Aj = — 1) and = if such a value of i 
does not exist. Re-combining the basis elements, we also can assume that the function / is 
either identically equal to zero or linearly independent with the functions , g^. There 
are equivalences between different subalgebras within the above classes induced by shifts and 
scalings of t. Thus, the argument t can be replaced in g and / by e^^i + eoi where eoi £1 £ In 
classes 4-6, a single nonzero parameter among cr, Aj, //j, Uj and k can be scaled to ±1. We also 
can set k = ±1 for classes 10 and 11 if k 7^ 0. An additional possibility for constraining signs of 
parameters is given by discrete symmetries of Eq. (jH). 

We omit the proof of the classification result because it is too cumbersome. 



4 Lie reduction and exact solutions 

As the vorticity equation (jj]) is an equation in three independent variables. Lie reduction using 
one- and two-dimensional subalgebras will lead to submodels of Eq. (jH) that are partial differ- 
ential equations in two independent variables and ordinary differential equations, respectively. 
Recall that it suffices to study only Eq. Q with 17 = 0. 

First we select inequivalent one-dimensional subalgebras of 5^, which are appropriate for 
Lie reduction. An optimal list of one-dimensional subalgebras of reads 

(P + aJi), {dt + aJi), (Ji + Z(5)), {Z{g)), 

where a G M and a G {—1, 0, 1} for the first and second cases, respectively, and g is an arbitrary 
smooth function of t, which does not vanish in the last case. See [Ij for more details. The algebra 
{Z{g)) cannot be used for classical Lie reduction as it does not permit making a reduction ansatz 
for ij:. This is why we only have to consider the first three subalgebras. 

For the subalgebra {V + aJi) the invariants of the associated one-parametric group are 
p = X — aint, q = n and v = tip. The reduction ansatz therefore is ■0 = t~^'v{p,q), which yields 

w + awp - VpWq + VgWp = 0, w:= y~2^pp + ~ (^) 

as a submodel of Though we have reduced the number of independent variables by one, it 
appears to be impossible to completely solve Eq. Solutions of Eq. ([8]) can be sought, for 
example, using numerical integration. We will return to a further reduction of Eq. ([8]) when we 
discuss reductions with respect to two-dimensional subalgebras. 

An ansatz constructed using the subalgebra {dt + aJi) \s = v{p, q), where p = X — at and 
q = This ansatz reduces Eq. @ to 

-{v + aq)qWp + {v + aq)pWq = Q, w := Y^-^Vpp + {{1 - q'^)vq)q. (9) 

Eq. ([9]) implies that w = F{v + aq), for F being an arbitrary smooth function of its argument. 
The definition of w then takes the form of a nonlinear Poisson equation. Li the particular case of 
F being a linear function, i.e. F = c{v + aq) + b, where b and c are real constants, this equation 
is a linear inhomogeneous PDE, 

-. ^ 2 Vpp + ((1 - q'^)vq)q = cv + caq + b. (10) 
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The general solution of Eq. (fTO|) is represented in the form v = Vg,h. + fp.i., where Up.i. is a 
particular solution of Eq. (jlOp and f g,h. is the general solution of the corresponding homogenous 
equation, which we denote by C. As Up.i. we can take the following solutions: 

b ca 

v = —q if c/-2,0, 

c c + 2 

6 , 2, 

v = Inl — (7 if c = 0, 

2 ' ^ ' 



b a 



V = - + -qlnll - q^\ if c = -2. 
2 3^' ^ ' 

In order to find solutions of C, we use the separation ansatz v = e*™^x(9)) where m is a 
complex constant and X is a complex- valued function of q. This ansatz reduces C to the ordinary 
differential equation 

2 



(TTl \ 
c + _ ^2 ) X = 0. 



The reduced equation is the Legendre equation whose general solution is expressed in terms 
of the associated Legendre functions and of first and second kind, where the parameter n 
is a root of the equation n(n + l) = — c, i.e. x = ^-fn*!?) with A and 5 being arbitrary 
complex constants. The values m, Xi ^ and B can be complex due to the linearity of C as the 
real part of any solution of C is also a solution of C. Combining solutions obtained for a fixed 
value of n and different values of m, we obtain, e.g. for c = — n(n + 1)7^ —2, 0, the solution 

Af 

v{p, «) = J] Re (A,e^'"^f iq) + P^e^'^^^'g™^ (g)^ ''''^ 
=1 



c + 2 



(If c ^ 0, the constant b can be neglected up to shifts of v, which are Lie symmetry transfor- 
mations for Eq. Q induced by shifts oi ip.) This gives rise to group-invariant solutions of the 
form 



N 

c + 2 

where we already take into account the effects of rotation by reverting to the moving coordinate 
frame using transformation 

For the above solutions to be globally defined on the sphere, the parameters involved should 
satisfy special conditions. Namely, we have that n G N, all rrij take only integer values from — n 
to n and the constants Bj are zero. In other words, the globally defined solutions of the kind 
considered are expressed in terms of spherical harmonics of the same degree n, 

^-^ n(n + 1) — 2 

where n € N and uj = \ — {a — Q)t, Am and Sm are real constants, is an associated Leg- 
endre polynomial, P^{^) = (— 1)™(1 — ^J?)"^/"^ {dJ^ / d^"^)Pn{^j) and P„ is the ordinary Legendre 
polynomial of degree n. The solutions of this form can be interpreted as traveling waves in 
East-West direction superimposed to a West-East mean flow that only depends on the latitude. 
Prom the physical point of view, this class of solutions is reasonable as it neither violates the 
cyclic boundaries on the sphere nor has singular points (as would be the case for non-integer 
values of m and n). Moreover, it describes a recurrent configuration for a barotropic atmosphere 
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at mid-latitudes and it can be found qualitatively on almost every weather chart for the height 
of the 500 Hectopascal pressure surface. 

In the more specific case when we require the mean flow to vanish, the relation a = 
must hold for the constants a and n. In this case, the phase velocity of the waves is given by 

Cphase — 7 ; TTj 

n[n + Ij 

which is the well-known phase velocity of Rossby-Haurwitz waves. This is probably the most 
famous solution of the barotropic vorticity equation on the sphere. Evaluating the formula ()lip 
for a few nonzero values of A^^s and a lower value of n already gives a quite nontrivial and 
physically reasonable form of the stream function tp, cf. Fig. [TJ 




Figure 1. Plot of the generalized Rossby-Haurwitz wave solutions for n = 20 at t = 0, where A3 = 100, 
53 = 0, y4g = 150, (Jg = 1.5, A\z = 200, &\z = 3.4 and Aig = 250, 5i8 = 0.9. The remaining Ai, Si are identically 
zero. The solution evolves against the direction of rotation of the Earth in westward direction. 

The reduction associated with the subalgebra (j7i + 2(g)) is of particular interest for the 
present paper letter. The reduction ansatz is -0 = v{p,q) + g{t)X, where p = t, and q = fi. 
The reduced form of Eq. 1^ is Wp + gwq = 0, w := ((1 — q^)vq)q. This equation is completely 
integrable by quadratures, which leads to the family of exact solutions 

i) = g{t)\ + f{t) + h{t)aiciaii\i^i + j I ^^^^^^ ^ ^ _ J g(j.^^^ (12) 

of Eq. dH), where f{t), h{t) and w{'y) are arbitrary functions of their arguments. 

Let us now briefly discuss two-dimensional Lie reductions of Eq. Q. An optimal list of 
two-dimensional subalgebras is 

{V + bJi,dt), {V,Ji + Z{at-^)), {V + aJi,Z{\t\')), 

{dt,Ji+Z{c)), {dt + cJi,Z{e'')), {Ji+Z{g),Z{g)), {Z{g'), Z{g^)) , 

where a,6 € M, c € { — 1,0,1}; c € { — 1,0,1} if c = 0; g, g, g^ and g'^ are arbitrary smooth 
functions of t such that g^ and g^ are linearly independent and g does not vanish. We first 
note that Lie reductions using the third, the fifth, the sixth and the seventh algebra are not 
possible as they include subalgebras of the algebra {Z{g)), which does not allow setting up a 
reduction ansatz. Then, the second and the fourth algebra include the subalgebra of the form 
{J'i + Z{g)) with the values g = at^^ and g = c, respectively. As we have seen above, the reduced 
equation stemming from the algebra {Ji + Z{g)) can be integrated completely by quadratures 
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for arbitrary g. Therefore, no further Lie reduction with respect to two-dimensional subalgebras 
containing (j7i + 2(g)) must be carried out. It thus remains to investigate the Lie reduction 
corresponding to the first subalgebra {V + aJi, dt). An ansatz ip = e~^^^v{fi) constructed using 
this subalgebra coincides, up to notation, with a particular case of the ansatz which is derived in 
the course of the reduction with respect to the subalgebra [dt + aJi) (with a = 0) and subsequent 
separation of variables in Eq. ([9]) with linear F. Summing up, two-dimensional Lie reductions 
of Eq. dH) give no new closed-form solutions in comparison with one- dimensional reductions. 

We also try to apply other tools for finding exact solutions to Eq. such as partially 
invariant reductions and hidden symmetries. 

In fact, the algorithm of partially invariant reduction cannot be directly used for Eq. (jH) 
since it is a single equation in a single unknown function, which is the stream function. At 
the same time, usually we can replace a single differential equation by an equivalent system of 
differential equations with greater number of unknown functions. There are at least two quite 
natural counterparts of the vorticity equation among systems of differential equations, namely, 
the corresponding Euler equations ([TJ for the flow velocity and the pressure and the vorticity 
equation itself, where the vorticity is assumed as one more unknown function and its definition 
is treated as an equation relating C, and if). We consider the last system, which we will simply 
call system dH). It is convenient that the maximal Lie invariance algebra of system is 
isomorphic to the algebra . More precisely, every vector field from is a prolongation of 
an operator from to the vorticity C, in view of its expression via the stream function ^. This 
is why the lists of subalgebras presented in Section [3] can be used for finding partially invariant 
solutions of system ([U . It suffices to consider only subalgebras whose intersections with the ideal 
{Z{g)) are equal to (^(1)). We have tested two partially invariant reductions, with respect to the 
subalgebras (j7i,^(l)) and {V + aj'i, dt, 2(1)) . After cumbersome computation we have derived 
that all the corresponding partially invariant solutions are in fact Lie invariant. Nevertheless, 
we plan to continue the study of partially invariant solutions for systems related to Eq. (jj]) since 
this may possibly result in interesting solutions for Eq. 

Hidden symmetries are defined as Lie symmetries for submodels [1], and only purely hidden 
symmetries, which are not induces by Lie symmetries of the initial model, are interesting. Un- 
fortunately, all Lie symmetries of both the valuable reduced equations ([8]) and Q are induced 
by symmetries of Eq. (jH). At the same time, the homogeneous linear counterpart C for the once 
integrated equation ([9]) admits, in addition to the induced Lie symmetry dp, the Lie symmetries 
vdy and h{p, q)d.u which arise due to the linearity of C and are purely hidden for Eq. @. Here 
the function h = h{p, q) traverses the set of solutions of the equation C The above symmetry 
allows us to separate variables and to linearly combine solutions. This is the fact which explains 
from the symmetry point of view why Lie reduction with respect to the class of subalgebras 
(dt + aJi) gives the widest and most interesting family of solutions. 

Let us now compare the reductions from |3] discussed above with those carried out in [8], 
where Lie reductions were computed for two particular cases of two-dimensional subalgebras. It 
turns out that these cases are equivalent to each other and, moreover, can be derived as special 
instances of the exact solution (|12p constructed in the course of Lie reduction with respect to 
the one-dimensional subalgebra {J'l -\- 2(g)). 

Explicitly, the first case is the Lie reduction with respect to the algebra similar to (P, J7i). 
The corresponding solution 

2i'Rq u t t 
where H'(6) = F(9), is included as a particular instance in the exact solution (|12p by setting 

g(t)=0, f(t) = ^, m = -% w(^) = 0. 
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Note that the first term in the solution tp^^^ is not essential as it accounts for the rotation of 
the reference frame which can be taken into account for the solution (jl2p in the nonrotating 
reference frame upon applying the transformation ([6]). The second term in tp^^'^ arises due to 
invoking the stream function only for the perturbation velocities and can be removed upon using 
the transformation ([5]). 

The second closed-form solution was constructed in [8] in the course of reduction using the 
algebra similar to (DjJ's) and it reads 

cos 9 H{6) Ci C2, I 



j,m = 1^ _ flSZl + ill + hl^ In 
^ 2uRo V t t 



1 + K 



where k = sin6'cos(A-|-t/(2i?o))- However, this solution is related to the above solution ■p^'^'^ and 
hence to the solution (jl2p via a point transformation. Namely, from the adjoint action 



Ad(e^^2 ) = _ J3 sin e + Ji cos e, 

it is clear that for the value e = — vr/2, the operator Ji is mapped to J^. The same transformation 
then maps the solution ()12p to the solution including the particular form ^/^t^l. Explicitly, the 
finite transformation component of J2 for fi is fl = ficose + y^l — /^^ cos A sin e, which for 
e = —tt/2 gives jl = — /U^ cos A. The transformed form of ()12p for g(t) = then takes the 

form 



ip = f{t) — /i(t)arctanh - cos A) + W{y^l — ^i^ cos A), 

where W = W{y/l — ^2 cos A) is an arbitrary function of its argument. The solution 'ijP'^ is the 
particular case of the above solution where 

f{t) = ^, h{t) = -^, W = 0. 

Again, the first two terms in -(/jl^l are inessential in that they can be added upon transforming 
to the rotating coordinate system and redefining the stream function. 

In short, the two exact solutions 'ip^'^'^ and V"'^^ and all the special forms therefrom studied in [S] 
can be united upon investigating the solution (|12p stemming from the algebra {J'i + Z{g)), which 
belongs to the list of inequivalent one-dimensional subalgebras of the maximal Lie invariance 
algebra ([7]). 



5 Concluding remarks 

The present letter is devoted to the relation of the two independent studies [3] and [S] on a 
particular form of the two-dimensional incompressible nonlinear Euler equations in a rotating 
reference frame in spherical geometry. The main intention of this work is to point out the impor- 
tance of systematically following the algorithm of Lie reduction in order to arrive at inequivalent 
exact solutions of a system of differential equations. The key aspects of this algorithms are the 
following: 

(i) Identify essential and inessential parameters in the system of differential equations at 
hand. This point is systematically addressed by solving the group classification problem, 
see e.g. \12\ I15j. In the case of the vorticity equation in a rotating reference frame ([ID the 
solution of the group classification problem gives that O can be set to zero upon using the 
simple point transformation Q. 

(ii) Determine an optimal list of inequivalent subalgebras of the maximal Lie invariance algebra 
of the system of differential equations of interest which are suitable for Lie reduction. This 
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task usually requires to classify low-dimensional inequivalent subalgebras under the adjoint 
action of a Lie group on its Lie algebra, see e.g. [21 [TTl |12] for the relevant techniques 
involved in this step. It is of utmost importance to do this classification, as the result will 
give information on what Lie reductions should be carried out for the particular system of 
differential equations. In the case of the vorticity equation (jH) a suitable optimal system 
was found in [1]. As no such optimal system was established in [8], the Lie reductions that 
were carried out are not guaranteed to give inequivalent solutions. Indeed, as is shown in 
Section H] the two-dimensional reductions studied in [8] are equivalent to each other with 
respect to Lie symmetry transformations of Eq. ([!]). 

(iii) If applicable, study the physical importance of the solutions obtained. The question of 
physical meaning usually cannot be answered solely in the framework of group analysis. 
However, what symmetries can add to this investigation is that any solution obtained in 
the course of Lie reduction is in fact only a representative of a parameterized family of 
similar solutions, which can be obtained upon acting with symmetry transformations on 
the derived solution. In some cases, an unphysical exact solution is simply the consequence 
of choosing an inappropriate ansatz for Lie reduction. 

Depending on the complexity of the maximal Lie invariance algebra of the system of differ- 
ential equations under consideration, the exhaustive solution of the Lie reduction problem can 
be a tedious task. It might therefore be tempting to compute only the Lie symmetries of a 
differential equation (which is an algorithmic task that can be handled with the aid of computer 
algebra) and then to use certain obvious combinations of the infinitesimal generators in order to 
carry out the reduction. Combinations that arise quite often for equations in sciences are those 
leading to scale-invariant or wave solutions. On the other hand, even if such solutions might 
have an immediate physical interpretation, there is no guarantee that they are not equivalent 
to each other. A systematic approach to Lie reduction therefore cannot avoid the classification 
of inequivalent subalgebras. 

Another issue that arises in the course of Lie reduction and which can considerably simplify 
the solution of the Lie reduction problem is whether one is able to completely integrate the 
equations arising at a low stage of reduction. In this case, the reduction with respect to higher- 
dimensional algebras that include as a subalgebra the respective algebra which enabled the 
integration can be omitted. This holds true in the case of the two-dimensional Lie reductions 
carried out in [^, which all follow from a one-dimensional reduction of the vorticity equation 
that can already be integrated completely. 

Lie reduction is only one method to derive exact solutions for systems of differential equations. 
The notion of point symmetries can be extended beyond point character and point transforma- 
tions can be used in more sophisticated methods. Other techniques that can lead to valuable 
physical solutions are e.g. conditional symmetries, potential symmetries, nonclassical symme- 
tries or the method of differential constraints. In order to properly apply these techniques, 
however, an exhaustive understanding of the possible Lie reductions is of certain value. This 
holds in particular for the techniques that are extensions of the classical Lie method, such as 
the method of differentially invariant solutions, the direct (ansatz) method, which is equivalent 
to the reduction method using conditional symmetries, or partially invariant reductions with 
Lie symmetries. Note that the latter have been tested for the vorticity equation on the sphere 
in this letter and it was found that they reduce to the group-invariant solutions that we found 
above. 

From the physical point of view, deriving exact solutions of atmospheric models in spherical 
geometry is a task of immediate relevance as such solutions can serve as important benchmark 
tests for numerical general circulation models. As an example for such benchmarks, in [18] a 
suite of seven tests was proposed for the assessment of numerical models for the shallow-water 
equations in spherical coordinates. By now, virtually every new numerical model of the shallow- 
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water equations is tested using the suggested single benchmarks, ah of which ahow the evaluation 
of different relevant aspects of the shallow-water equations. By providing a host of exact and 
physically relevant solutions for the barotropic vorticity equation in spherical geometry, with 
this letter we aim to extend the list of known solutions of models of geophysical fluid dynamics 
on the rotating sphere. We should like to stress in particular that using the method of Lie 
reduction, several exact solutions can be derived that extend or generalize known solutions of 
the vorticity equation, such as the well-known Rossby-Haurwitz wave solution. It is therefore 
hoped that the general classes of exact solutions found will be of further relevance in the field 
of numerical geophysical fluid dynamics. 
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